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This  is  a  final  report  on  AFOSR-78-3533  and  includes  three  attachments: 
Report  SM  80-14,  "An  Improved  Model  for  Second  Phase  Particle  Interaction 
Analysis,"  by  R.R.  Batarags;  Report  SM  81-1,  "A  Simple  yet  Effective 
Substructuring  Procedure  for  Finite  Element  Analysis,"  by  M.  Kondo  and  G.B. 
Sinclair;  and  an  abstract  of  the  paper,  "Computation  of  Slow,  Stable  Crack 
Growth  in  a  Fracture  Test  Specimen  using  the  J-Integral,"  by  J.L.  Swedlow, 
F.B.  Prinz,  Z.-C.  Lai,  and  C.W.  Cho. 

The  Research  Grant  covered  the  period  1  November  1977  to  31  October  1981 
and  was  entitled,  "Analysis  of  Critical  Microscopic  Flew  Behavior."  The 
research  objectives,  in  general  terms,  were  to  simulate  events  preceding 
fracture  in  selected  alloy  systems  via  computational  plasticity;  the  aim  was 
to  quantify  these  events  beyond  what  is  possible  using  standard  metallurgical 
tools  (e.g.,  microscopy,  optical  and  electronic).  In  particular,  it  was  held 
that  interactions  between  second  phase  particles  would  significantly  alter  the 
size  scale  of  damage  in  high  strength  alloys.  [A  detailed  recapitulation  of 
this  hypothesis  appears  cm  pp.  1-2  in  the  Batarags  report.] 

At  the  outset,  analyses  were  performed  to  this  end.  While  the  results 
tended  to  confirm  the  hypothesis,  subsequent  review  of  the  models  used  (i.e. , 
the  finite  element  maps)  shewed  that  the  models  themselves  were  deficient  in 
that  their  design  tended  to  force  the  desired  result.  [A  critique  of  these 
early  results,  together  with  an  extensive  sampling  of  the  corresponding  data, 
appears  on  pp.  3-8  of  the  Batarags  report.] 

Next,  Batarags  designed  and  throughly  tested  a  new  model  which  met  the 
criteria  that  we  had  begun  to  recognize  were  crucial  for  the  successful 

resolution  of  this  problem.  It  is  believed  that  Batarags  was  successful  in 


this  respect,  and  his  report  is  concentrated  on  those  findings.  In  that 
respect  and  in  others  referred  to  above,  this  report  stands  as  the  central 
technical  docusent  for  this  Research  Grant. 

At  this  stage  the  hunan  element  intrudes.  Mr.  fiatarags  finished  his 
studies  for  the  Master  of  Engineering  degree  and  left  the  University.  Also, 
the  Principal  Investigator,  Professor  Swedlow,  had  been  found  to  have  a  brain 
tunor  and  was  scheduled  for  appropriate  surgery  in  the  Spring  of  1981.  The 
upshot  was  that  the  personnel  who  had  been  engaged  in  the  mainstream  of  work 
on  this  Research  Grant  were  unable  or  unavailable  to  continue  it.  In  effect, 
no  reportable  work  on  the  matter  of  particle  interaction  was  performed  beyond 
the  date  of  the  Batarags  report. 

It  happened,  however,  that  a  new  student,  Mr.  Lai,  became  available  late  in 
the  Spring  of  1981,  and  the  PI  had  recovered  sufficiently  to  resume  his 
duties.  After  a  modest  learning  period  during  which  the  new  student 
assimilated  the  prior  work,  an  effort  to  continue  was  made.  In  the  relatively 
short  time  then  available  prior  to  the  end  of  the  granting  period,  no  further 
reportable  work  was  concluded  on  this  aspect  of  the  research.  Since  that 
time,  however,  additional  analyses  were  concluded  and  a  report  is  now  in  draft 
form.  Its  abstract  is  attached. 

However,  another  line  of  inquiry  was  in  progress.  Professor  Sinclair,  the 
Associate  Investigator,  had  been  examining  the  issue  of  debonding  between  a 
particle  and  the  surrounding  matrix  material.  Crucial  to  this  work  is  the 
question  of  what  condition  governs  the  event.  His  approach  was  to  hypothesize 
a  local  criterion  and  then  to  determine  whether  the  results  were  in  some  sense 
acceptable.  At  the  outset,  for  example,  he  tested  the  notion  of  using  a  local 
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critical  stress.  If,  at  sane  part  of  the  boundary,  debonding  began,  it  became 
evident  that  the  initial  debond  would  form  a  crack;  this  configuration, 
however,  is  such  that  debonding  must  continue  owing  to  the  large  local 
stresses  induced  at  the  crack's  tip.  Clearly  this  forces  the  debonding 
process  to  became  unstable,  and  particle-particle  interaction  would  be  greatly 
enhanced.  Since  this  sort  of  instability  is  not  commonly  observed,  it  was 
concluded  that  a  local  stress  criterion  for  debonding  is  not  wholly 
appropriate,  and  a  local  energy  release  rate  GA  was  put  forward  as  the 
quantity  governing  debonding.  This  condition  was  not  fully  tested  prior  to 
the  end  of  the  granting  period. 

To  follow  his  line  of  inquiry,  Professor  Sinclair  concluded  that  a  fine 
degree  of  resolution  would  be  required.  Normally,  this  would  equate  to  very 
large,  very  detailed  models.  Such  finite  element  maps  are  infeasible  to  work 
with  on  the  University's  computers,  and  he  therefore  developed  a  technique  for 
substructuring  described  in  the  attached  report.  This  is  a  useful  and 
effective  technique  for  the  purpose  intended  but  was  not  followed  up  with 
further  studies  of  the  debonding  problem. 

The  status  of  the  research  performed  under  this  grant,  therefore,  is  that: 

1.  The  research  questions  have  been  significantly  sharpened; 

2.  the  requisite  (computational)  tools  have  been  developed  and 
dooueented;  and 

3*  there  has  been  uisponsored  follow-on  work  now  in  the  process  of 
docwentation 

It  la  regrettable  that  the  research  objectives  have  not  been  met  during  the 
granting  period  but,  in  retrospect,  the  research  objectives  are  now  seen  as 
overly  ambitious  and  will  probably  take  a  number  of  years  further  to  achieve. 
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Owing  to  the  status  of  this  work,  none  of  the  doc mentation  has  so  far  been 
submitted  for  publication.  The  report  in  preparation,  item  3  above,  will  in 
all  likelihood  be  published. 

In  the  same  vein,  there  have  been  no  coupling  activities  nor  inventions  on 
this  Research  Grant. 
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I.  INTRODUCTION 


The  rol*!  of  microstructure  and  how  its  behavior  under  loading  influences 
the  ductile  fracture  process  are  of  significant  interest  to  metallurgists  as 
well  as  the  mechanics  comnunlty.  We  have  addressed  the  issue,  centering  our 
efforts  on  computer  simulation  of  mlcrostructural  events  preceding  fracture  - 
primarily  dealing  with  the  effect  of  second  phase  particles.  To  date  our  efforts 
have  met  with  limited  success.  However,  we  believe  that  the  modeling  of  the 
microstructure  under  consideration  which  comprises  the  bulk  of  this  report, 
has  been  improved  and  should  allow  for  further  work  to  proceed  at  a  more  rapid 
rate. 

In  order  to  provide  a  framework  for  our  efforts,  as  well  as  to  impart  some 
sense  to  the  sections  to  follow,  a  recapitulation  of  how  we  envisage  the  ductile 
fracture  process  to  proceed  Is  presented  first. 

1.  A  dispersion  of  second  phase  particles,  whose  material  properties  differ 
from  the  bulk  of  the  material  (here  termed  the  matrix  and  which  is  considered 
both  homogeneous  and  Isotropic)  will,  under  the  action  of  low  level  loading, 
cause  a  local  stress  concentration  to  appear. 

2.  As  the  level  of  load  Increases  and  small  scale  (i.e.  contained) 
plastic  yielding  occurs,  the  elastic  stress  concentration  tends  to  be  suppressed, 
whereas  the  local  strains  (measured  In  terms  of  e.g.  the  octatedral  strain)  will 
tend  to  increase. 

2a.  The  foregoing  sequence  can  be  observed  in  the  vicinity  of  a  single 
second  phase  particle,  since  the  events  are  local  to  the  particle-matrix  interface 
not  long  range. 
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3.  Depending  on  variables  such  as  relative  material  properties,  particle- 
particle  spacing,  load  direction  vis-a-vis  the  particle-particle  axis  and  lateral 
constraint,  a  level  of  excitation  will  be  attained  at  which  the  size  of  the  plastic 
zones  are  no  longer  confined  to  and  determined  by  an  individual  particle,  but 

will  rather  connect  with  the  plastic  zones  surrounding  a  neighbor  particle.  In 
other  words,  a  change  of  size  scale  can  occur  such  that  the  extent  and  shape  of 
the  plastic  zone  is  no  longer  governed  by  a  characteristic  dimension  of  a  single 
particle,  but  is  instead  governed  by  a  dimension  relating  to  the  proximity  of 
neighboring  particles. 

4.  Once  plastic  zones  from  neighboring  particles  have  interacted,  further 
excitation  will  tend  to  concentrate  within  such  highly  strained  regions.  The 
interaction  of  second  phase  particles  creates  a  pathway  for  macroscopic  ductile 
fracture.  Further  material  damage,  such  as  the  retail urgi cal ly  observed  void 
growth  and  coalescense  phenomena,  will  occur  in  the  most  severely  strained 
regions,  leading  to  development  of  small  fissures  initially,  and  gross  section 
fracture  ultimately. 

Our  efforts  have  been  concentrated  on  the  effects  of  interacting  second 
pnase  particles,  and  on  the  development  and  extent  of  strain  localization 
We  have  attempted  to  demonstrate  that  the  mere  proximity  of  neighboring  particles 
can  lead  to  interaction  between  the  attendant  plastic  zones,  at  quite  low  far 
field  levels  of  excitation.  Furthermore  and  most  importantly,  we  have  attempted 
to  show  that  the  presence  of  multiple  second  phase  particles  will  accelerate 
the  growth  of  regions  of  high  strain. 

Unfortunately,  we  have  of  late  become  sidetracked  from  the  main  issue  due 
to  modelling  problems,  which  are  discussed  in  the  next  section.  Consequently, 


the  microstructure  has  been  developed,  so  that  progress  can  be  made  on  the 


Initial  objective. 

II.  DISCUSSION  AND  EVALUATION  OF  PREVIOUS  WORK 


The  work  discussed  in  this  section  has  been  documented  earlier  [1,2].  The 
objective  here  is  not  to  detail  or  reiterate  the  findings  discussed  previously, 
rather  to  highlight  some  of  these  findings  for  continuity.  Additionally  to  be 
discussed  are  the  problems  with  those  findings  and  substantiating  evidence. 

The  initial  work  [1]  addressed  the  question  of  second  phase  particle 
interaction,  under  the  assumption  that  the  particle-matrix  interface  remained 
"coherent."  An  attempt  to  estimate  the  effect  of  particle  destruction  was 
made,  by  allowing  one  particle  to  open  up  at  its  center  (the  "burst"  model), 
although  such  a  procedure  was  acknowledged  to  be  not  entirely  realistic.  The 
latter  work,  [2],  concerns  itself  with  the  issue  of  a  matrix-particle  decohesion, 
and  the  effect  of  such  an  event  on  the  interaction  of  particles. 

The  particular  two  phase  system  chosen  for  these  analyses  is  a 
Titanium  alloy  (T1-6AL-4V),  consisting  of  cylindrical  second. phase  particles  (a) 
<n  an  otherwise  homogeneous  matrix  (s)  phase.  Elastically  then  phase  is  stiffer 
than  the  8,  however,  the  a  possesses  a  lower  yield  stress.  Furthermore,  the 
phase  exhibits  an  extremely  flat  stress-strain  curve  plastically.  That  is, a 
work  hardens  to  a  much  lesser  extent  than  the  e  phase.  The  stress-strain 
curves  for  both  the  a  and  8  phase  are  presented  in  Figs.  1_  and  respectively. 


The  fact  that  the  second  phase  particles  are  regarded  as  cylindrical  allows 
one  conceptually  to  section  the  material  in  such  a  fashion  as  to  achieve  a  two 
dimensional  problem  containing  circular  second  phase  particles.  Since  interaction 
is  the  issue,  a  minimum  of  two  particles  is  necessary.  In  these  studies,  the 
center  to  center  distance  was  held  fixed  at  11  1/2  diameters.  Because  this 
alloy  system,  heat-treated  differently,  can  produce  a  wide  variation  in  particle 
spacing,  11  1/2  diameters  is  very  conservative. 

Both  analyses  utilized  the  identical  finite  element  model  for  the  microstru.ture, 
which  is  presented  in  Figures  3-5.  Figure  3  depicts  the  complete  model  allowing 
for  two  second  phase  particles  and  three  of  the  five  distinct  orientations  of  the 
particle-particle  axis  with  respect  to  the  loading  axis,  assumed  vertical.  Figure 
4  indicates  the  regions  encompassing  both  second  phase  particles,  and  their 
connectivity.  An  enlargement  of  one  particle,  and  its  immediate  surroundings 
is  seen  in  Figure  5. 

Excitation  of  the  model  (Figure  3)  is  effected  by  displacement  boundary 
conditions  -  tensile  in  the  vertical  direction.  Loading  in  the  lateral  direction 
is  not  as  easily  resolved.  At  least  three  distinct  lateral  loading  conditions 

(constraints)  can  be  envisaged: 

\ 

1)  no  motion  normal  to  the  direction  of  tensile  loading  permitted 
(Full  Constraint) 

2)  no  lateral  constraint  whatsoever  (Nil  Constraint) 

3)  lateral  motion  associated  with  Poisson  contraction  (Poisson  Constraint) 

A  fuller  discussion  can  be  found  in  SM77-7,  pages  7-9  (ref.  1). 

See  also  the  generalized  discussion  by  Swedlow  and  Smelser  (ref.  3). 


DISCUSSION  OF  SM77-7 


This  analysis  attempted  to  establish  the  phenomena  of  second  phase  particle 
Interaction  computationally,  at  far  field  strain  levels  smaller  than  those  typically 
viewed  by  metallurgists  studying  the  identical  problem.  The  intent  was  to  demonstrate 
that  the  presence  of  neighboring  particles  would  alter  the  strain  fields,  even  at 
such  low  levels  of  straining.  If  this  were  the  case  then  the  existence  of  second 
phase  particles  provides  a  precursor  to  large  scale  material  damage.  In  effect  - 
from  the  relative  material  properties  of  matrix  and  particles  and  their  separation  - 
the  load  under  which  thejnaterial  will  fail  ductidely  could  be  inferred. 

Representative  results  are  given  by  Figures  6-9.  The  first  two  figures 
(numbers  6  and  7)  show  octahedral  strain  values  around  one  particle.  The  data 
are  normalized  with  respect  to  octahedral  strain  at  the  proportional  limit,  a 
value  of  0.012257.  The  radial  location  is  in  the  matrix  material,  very  close  to 
the  particle-matrix  interface.  The  four  curves  indicate  varying  far- field  levels 
of  strain,  with  the  lowest  curve  being  an  elastic  result.  Figures  8  and  9 
depict  the  regions  in  which  the  local  octahedral  strains  exceed  far  field 
values  by  5%,  (i.e.  1.05  times  far  field  strain). 

At  the  time  the  results  were  obtained,  they  were  viewed  as  disappointing, 
although  not  inaccurate.  The  elastic  curves  on  Figures  6  and  7  indicate  no 
interaction,  as  should  be  expected  at  such  a  large  particle  spacing.  Furthermore, 
the  curves  are  smooth  and  reminiscent  of  elastic  results  for  say  a  hole  in  a 
plate.  From  Figures  8  and  9,  which  represent  the  least  and  the  most  interaction 
seen,  respectively,  again  benign  conclusions  were  drawn.  Consequently,  work 
proceeded  with  an  attempt  to  deal  with  matrix-particle  decohesion,  as  detailed 
in  the  next  section. 
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DISCUSSION  OF  SM79-8 

Using  the  identical  finite  element  model,  an  attempt  was  made  to  evaluate 
the  effects  of  a  local  matrix-particle  decohesion  on  the  ability  of  second 
phase  particles  to  interact.  The  so-called  burst  model  indicated  that  the 
regions  of  high  strain  would  exhibit  a  propensity  to  join  at  much  lower  levels 
of  far  field  loading.  Consequently,  the  assumption  was  that  a  decohesion  would 
also  accelerate  the  interactive  effects;  but  would  the  level  of  load  at  which 
such  a  decohesion  occurs  be  a  major  variable? 

Representative  results  are  given  by  Figures  10-16.  Figures  10  and  11 
show  the  variation,  with  angular  position  around  the  particle  at  which  the 
decohesion  was  enforced,  of  cR  Figures  12  and  13  -  the  variation  of  the 
normalized  octahedral  shear  stress,  both  sets  for  two  orientations.  Figures 
14-16  Indicate  regions  exhibiting  a  high  level  of  octahedral  strain  in  relation 
to  the  far  field  values. 

A  few  words  of  explanation  are  in  order  at  this  point:  Since  we  were 
and  still  are  lacking  a  criterion  for  decohesion,  the  decohesion  was  enforced 
at  different  levels  of  far  field  strain.  "Coherent"  model  implies  no  decohesion 
"Initial  debond"  Implies  particle-matrix  surface  is  separated  but  not  opened 
prior  to  application  of  any  load;  "debond  on  10"  -  small  local  plastic  yielding 
present  prior  to  debonding;  "debond  on  20"  -  plastic  flow  throughout  the  model 
prior  to  debonding;  "debond  on  15"  -  intermediate  load  level  to  "debond  on  10 
and  20."  Additionally,  the  extent  of  the  debond  is  approximately  1/12  of  the 
circumference  of  the  second  phase  particle. 

The  conclusions. drawa  from  the  analysis  were  that  the  orientation  of  the 
particle-particle  axis  with  respect  to  the  load  axis  was  a  very  important 
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varlable,  lateral  constraint  less  so.  Futhermore,  the  more  plastic  flow  present 
at  the  time  of  debond,  the  more  pronounced  the  tendency  for  the  regions  of  high 
octahedral  strain  to  propogate.  Finally,  that  the  level  of  interaction  seen  for 
any  analysis  with  a  decohesion  introduced  was  substantially  greater  than  the  level 
of  interaction  seen  in  either  the  corresponding  "cohesive"  or  "burst*'  model  case. 

DISCUSSION  OF  THE  MODEL  AND  ITS  LIMITATIONS 


The  finite  element  model  discussed  so  far,  shown  in  Figures  3-5,  has  severe 
limitations  for  the  purposes  intended,  although  this  is  not  obvious  from  the 
results  presented  and  discussed  previously.  Elastic  analyses,  such  as  can  be 
inferred  from  the  lowest  curves  on  Figures  6  and  7  are  not  indicative  of  the 
quality  of  results  obtainable  from  the  model  for  a  plastic  analysis.  That, 
however,  is  essentially  due  to  the  nature  of  plastic  flow  Itself.  Once  a 
region  has  yielded,  further  straining  tends  to  accrue  there,  as  opposed  to 
necessarily  yielding  the  surrounding  regions.  In  a  finite  element  sense,  an 
element  is  considered  as  yielded  if  the  yield  criterion  is  exceeded  at  the 
centroid  of  a  constant  strain  triangle  i.e.  a  binary  decision.  Consequently, 
if  a  large  element  extends  into  regions  of  both  high  and  low  stress,  a  good 
possibility  exists  of  its  not  yielding  even  if  nearly  50%  of  the  region  is 
experiencing  stress  levels  high  enough  to  initiate  plastic  flow.  Such  a  result 
may  appear  in  Figures  8  and  9  since  there  is  no  apparent  reason  for  the  particular 
distribution  of  elements  not  exhibiting  extensive  strains,  save  for  their  having 
yielded  later  than  their  neighboring  smaller  elements. 

The  bipolar  coordinate  scheme  used  to  facilitate  the  modeling  introduces 
several  other  geometric  flaws.  The  element  size  gradations  that  exist  are 
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prejudlclal  since  along  the  particle-particle  axis  the  elements  remain  rather 
small  between  the  particles,  but  not  as  much  in  other  directions  emanating  from 
the  particles.  In  effect,  the  plastic  flow  has  been  forced  to  occur  between 
the  particles.  Similarly,  there  Is  also  an  angular  bias  built  in  around  each 
particle  -  that  being  that  the  smallest  elements  around  the  particle  occur  directly 
on  the  particle-particle  axis  closest  to  the  neighbor  particle.  In  short,  the 
region  most  likely  to  yield  first  occurs  between  both  particles,  and  the  path  of 
least  resistance  to  plastic  flow  lies  directly  towards  the  neighbor  particle. 

The  existence  of  element  size  gradations  In  the  direction  of  the  neighbor 
particle  Is  obvious  from  Figure  4,  as  Is  the  angular  size  distribution.  Whether 
such  gradations  are  fatal  to  the  analyses  performed  is  definitively  answered  by 
Figures  17-19.  Th-se  three  figures  are  directly  analogous  to  Figures  14-16  with 
the  exception  that  the  results  shown  In  Figures  17-19  are  taken  from  analyses 
where  matrix  material  properties  were  substituted  In  the  region  normally  associated 
with  the  neighbor  particle.  That  Is,  Figures  17-19  are  for  a  one-particle  model. 

It  should  be  mentioned  that  even  though  these  results  show  regions  over  which  the 
normalized  octahedral  shear  strain  exceeds  far  field  values  by  5%  (Figures  14,  15, 
17,  and  18)  and  10%  (Figures  16  and  19)  the  entire  region  has  yielded.  Quite 
clearly  then,  if  Figures  14  and  17,  15  and  18,  16  and  19  are  viewed  as  pairs, 
the  findings  -  with  some  small  local  perturbations  -  are  identical!  Consequently, 
the  conclusion  must  necessarily  be  that  the  analyses  have  not  been  overly 
sensitive  to  the  presence  of  a  second,  second-phase  particle  with  second  phase 
material  properties,  but  rather  the  results  appear  dependent  and  quite  sensitive 
to  the  finite  element  model  itself.  Hence  the  model  must  be  rejected! 
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NEW  MODEL  AND  INDICATIVE  RESULTS 


Having  shown  the  bipolar  modeling  to  be  unacceptable  In  the  previous 
section,  we  embarked  on  devising  a  new  model.  A  number  of  considerations 
serve  both  as  guides  and  constraints.  The  requirements  and  limitations  we 
must  address  are: 

1)  A  clearly  defined  1:1  check  case  comparing  single  particle  results 
to  multiple  particle  results  must  be  available 

a)  In  a  model  sense  this  Is  easily  conceivable  and  must  be  easily 
obtainable 

b)  on  an  exltatlon  level  the  question  that  arises  and  still  needs 
a  definite  answer  Is  -  Is  the  comparison  made  on  equal  stress 
levels  or  equal  far  field  strain  levels? 

2)  Minimization,  If  not  elimination  of  element  size  gradations.  There 
are  two  distinct  gradations,  l.e.,  angular  and  radial,  both  of  which 
need  to  be  kept  In  mind. 

3)  The  ability  to  vary  orientation  of  the  particle-particle  connection 
axis  vis-a-vis  the  load  axis  must  exist. 

4)  The  ability  to  vary  L/D  (distance  between  particles  to  diameter  ratio) 
must  be  easily  attainable. 

5)  Points  1-4  to  be  met  without  excessive  computational  cost. 

The  scheme  settled  upon  attempts  to  consider  all  of  the  requirements  established, 
and  has  done  so  with  some  level  of  success. 

In  order  to  keep  the  model  size  within  limits  and  preserve  the  ability  to 
vary  L/D,  a  simple  coarse  to  fine  substructure  technique  Is  employed.  The  coarse 
or  matrix  model  Is  shown  in  Figure  20,  where  the  cell  along  which  displacements 
are  carried  over  to  the  fine  or  particle  model  is  indicated  by  the  heavy  point. 
Figures  21  and  22  Indicate  how  the  L/D  ratio  can  be  varied  from  two  to  four 
approximately.  The  L/D  Is  approximate  as  each  particle  is  modeled  as  a  square, 
whose  area  is  then  converted  into  a  circle  of  equivalent  area  yielding  a  diameter 
(D).  L  is  measured  from  center  to  center. 


A  single  particle  problem  is  analyzed  by  substituting  second  phase  material 
properties  into  only  the  center  square  region  darkened  in  Figure  22.  The  single 
particle  case  is  the  reference  case,  to  which  all  further  multiple  particles 
cases  must  be  compared.  The  comparison  is  imperative  since  the  question  being 
asked  is:  to  what  extent  are  the  high  strain  regions  altered  due  to  the  presence 
of  neighboring  particles? 

The  particle  model  is  presented  in  Figures  23  and  24  in  its  actual  size 
and  configuration,  and  an  enlargement,  respectively.  The  model  as  presented  can 
easily  offer  two  distinct  length  to  diameter  ratios  i.e.  L/D  *  2  and  L/D  *  4. 
Furthermore,  larger  length  to  diameter  ratios  can  be  designed  by  simply  decreasing 
the  size  of  the  particle  itself  on  the  matrix  model.  This  does,  however,  actually 
create  new  finite  element  models,  which  we  have  chosen  not  to  depict  at  this 
time.  In  fact,  we  have  dealt  solely  with  the  L/D  =  2  model,  trying  to  ensure 
that  it  performs  sufficiently  well  for  our  needs.  The  close  spacing  is  a  result 
of  a  change  of  perspective.  Since  the  previous  analyses  dealt  with  L/D~ll  1/2, 
and  were  somewhat  unsuccessful,  it  was  decided  that  a  more  reasonable  approach, 
especially  considering  the  difficulties  encountered  with  the  model  itself,  is 
to  analyze  a  situation  where  heavy  interaction  is  certain  to  occur.  Hence  the 
close  spacing. 

From  Figure  24,  it  is  obvious  that  angular  size  gradations  do  not  exist. 
However,  modest  radial  size  gradients  are  present.  Essentially  this  is  due  to 
computer  limitations.  We  recognize  and  acknowledge  that  the  model  is  a  compromise 
However,  we  also  believe  it  to  be  vastly  superior  to  the  previous  model. 

Variations  of  the  particle-particle  axis  with  respect  to  the  load  axis 
has  been  approached  from  a  slightly  different  viewpoint  than  previously.  The  now 
discarded  model  actually  consisted  of  five  different  models,  three  of  which  are 


given  by  Figure  3.  We  have  chosen  to  effect  the  rotation  of  load  axis  with 
respect  to  particle-particle  axis  through  the  input  boundary  conditions.  That 
is,  instead  of  imposing  uniaxial  tension  on  the  model  shown  in  Figure  20  vertically, 
we  can  Impose  boundary  conditions  to  be  those  of  say  45°  off  vertical  on  the 
model.  Through  transformation  equations  any  off  angle  (from  vertical)  boundary 
conditions  can  be  calculated. 

The  microstructure  we  are  modeling  contains  cylindrical  second  phase  particles, 
which  should  be  modeled  as  circles  in  two  dimensions.  Since  the  matrix  model. 

Figure  20,  models  the  second  phase  particles  as  squares,  some  justification  is 
in  order.  During  the  attempts  to  create  a  new  model,  we  had  designed  a  very 
similar  model  to  that  presented  here,  although  containing  somewhat  less  detail. 

Using  that  model  we  analyzed  the  problem  of  a  hole  in  an  elastic  plate  in  tension. 
The  problem  was  solved  in  two  distinct  ways  that  are: 

1)  direct  loading  of  the  fine  model,  and 

2)  loading  of  the  matrix  model,  and  subsequent  substructure  loading,  i.e., 
of  the  fine  model. 

The  Importance  Is  two-fold.  First  we  can  compare  results  to  well-known  analytic 
solutions  and  secondly  loading  of  the  matrix  model  and  subsequent  loading  of  the 
particle  model  Informs  us  of  the  quality  of  results  we  can  obtain  by  an  extremely 
crude  modeling  of  a  circular  particle  by  a  square.*  Both  analyses  showed  excellent 
agreement  with  theory.  Since  the  material  properties  of  the  second  phase  are  not 

*When  the  boundary  conditions  are  Input  on  the  matrix  model,  a  square  region 
similar  to  that  seen  on  the  center  of  Figure  20  possesses  essentially  nil 
elastic  properties.  The  displacements  are  then  carried  onto  the  particle  model 
where  a  circular  hole  is  vastly  better  modelled. 


nearly  as  dissimilar  form  the  matrix  as.  compared  to  nil  elastic  properties,  we 
expect  two  material  problems  to  show  the  same  kind  of  agreement.  In  short,  we 
feel  justified  in  our  approximation  of  a  circular  region,  on  the  matrix  model, 
to  be  a  square. 

DISCUSSION  OF  RESULTS  -  Testing  of  the  Element  Model 

We  have  run  a  rather  extensive  list  of  elastic  check  cases  which  essentially 
are  an  attempt  to  insure  that: 

1)  we  have  no  built-in  bias 

2)  we  have  sufficient  detail  on  the  particle  model 

3)  we  are  justified  and  introduce  no  difficulties  by  rotating  our  input 
to  simulate  rotation  of  the  particle-particle  axis  WRT  the  load  axis. 

We  also  have  preliminary  plastic  results  available  which  will  be  presented  as 

well. 

Although,  from  our  experience,  elastic  results  are  not  necessarily  conclusive 
indicators  of  the  quality  of  plastic  results,  it  is  nevertheless  important  to 
examine  elastic  results  in  some  detail.  Furthermore,  with  the  present  model,  we 
can  obtain  a  sense  of  the  quality  of  our  model,  since  we  can  compare  single 
to  multiple  particle  results.  Additionally,  we  expect  to  see  differences  in  the 
elastic  results  since  L/D  »  2.  Therefore,  a  comprehensive  set  of  elastic  results 
are  presented  in  Figures  25  and  26.  The  angular  variation  of  for  various 
locations  are  given  by  Figures  25-32,df  <j00/a  by  Figures  33-40, of  Tr0/a  by 
Figures  41-48.  Normalized  octahedral  shear  stress  angular  variations  are  depicted 
by  Figures  49-56.  The  radial  location  measured  from  the  center  of  the  particle 
starts  at  the  closest  element  ring  in  the  matrix  material  and  progresses  outwards 
for  each  of  the  sets  of  data  presented. 


The  results  presented  substantiate  our  belief  that  the  present  model  is 
a  vast  improvement  over  the  previous  model.  The  darkened  symbols  indicate  single 
particle  results  for  various  angles  of  far  field  tension  with  respect  to  particle- 
particle  axis.  Since  these  points  fall  essentially  onto  one  single  curve,  we 
can  safely  conclude  that  our  method  of  simulating  differing  particle-particle 
axis  orientations  WRT  the  loading  axis  through  input  boundary  conditions  introduces 
no  discernable  errors. 

Since  these  analyses  were  performed  for  a  particle  spacing  given  by  L/D=2 , 
we  expect  to  see  differences  in  the  stress  components  once  we  analyze  the  multiple 
particle  cases.  These  differences  do  appear  as  well.  In  general,  from  Figures 
49-56  we  can  say  that  the  stress  level  is  somewhat  higher  i.e.,  clearly  distinguishable 
from  the  single  particle  results.  The  scatter  of  the  data  is  somewhat  greater  for 
the  multiple  particle  results,  than  the  single,  but  not  objectionably  so. 

A  rather  interesting  event  is  noticeable  on  Figures  49-56  for  the  multiple 
particle  data.  The  nominal  average  value  for  Toct/Tiim1s  on  the  order  *28, 
maximum  “.32.  Now  as  we  proceed  radially  outwards  from  tne  particle,  the  trend 
that  appears  is  for  the  peak  values  of  tq  to  remain  roughly  constant, 

while  the  minimum  values  tend  to  increase.  Certainly  the  fact  that  the  curves 
tend  to  flatten  out  is  not  surprising  -  in  fact  for  single  particle  results  at 
a  large  enough  distance  we  expect  the  Toct/Ti im  curve  t0  be  Nevertheless, 

why  do  the  'valleys'  get  deeper  in  Figure  49,  rather  than  the  peaks  becoming 
higher  in  comparison  to  Figure  56? 

For  our  purposes  the  more  important  question  is  how  does  this  particular 
model  perform  when  plastic  yielding  occurs?  Preliminary  results  have  been  obtained 
and  are  presented  in  Figures  57-62.  These  figures  graphically  depict  regions 
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of  high  strain  (measured  by  the  octahedral  shear  strain)  for  both  single 
particle  and  multiple  particle  analyses.  The  figures  represent  data  obtained 
for  an  analysis  where  the  far  field  loading  is  presumed  tension  at  22.5°  from 
vertical  on  the  matrix  model,  at  an  accumulated  strain  level  of  “4 %.  There 
are  very  distinct  differences  between  single  and  multiple  particle  results. 

On  Figures  57  and  58  the  elements  exceeding  far  field  strain  levels  by  5%  are 
darkened  and  a  shift  Is  apparent  from  extensive  yield  in  the  load  line  direction 
(Figure  57)  to  preferred  yield  towards  the  neighbor  particle  (Figure  58).  In 
fact  on  the  basis  of  the  105%  characterization,  the  particles  have  already 
connected.  As  we  use  a  stricter  characterization  -  107%  for  Figures  59  and  60, 
and  110%  for  Figures  61  and  62,  of  course  the  high  strain  regions  become  pro¬ 
gressively  smaller.  However,  in  all  cases  the  effect  of  multiple  particles 
is  clearly  visible. 

One  note  of  caution  should  be  mentioned  concerning  the  results  given  on 
Figures  57-62.  The  results  are  not  precisely  symetrical  as  we  have  depicted 
the  data.  The  principal  reason  is  since  the  tension  axis  is  at  22.5°  there 
is  no  precise  symetry  of  element  distributions.  For  example,  had  we  chosen 
to  have  the  tension  axis  at  a  multiple  of  30°,  the  elements  would  be  distributed 
precisely  symetrically  about  the  tension  axis.  Clearly,  a  polygonal  representation 
for  a  circle  has  a  finite  number  of  such  symetry  lines.  Nevertheless,  had  the 
load  been  applied  in  such  a  symetrical  direction,  we  would  expect  symetrical 
results.  In  addition,  the  use  of  a  5%  exceedance  of  the  strain  criteria  is 
rather  arbitrary  and  does  not  account  for  the  minor  numerical  variations  which 
occur  in  any  analysis  of  the  type  presented. 
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An  Interesting  observation  can  be  made  by  viewing  Figures  62,  60  and  58, 

In  that  order.  Although,  we  are  then  viewing  a  progression  of  less  severe 
criteria  for  the  determination  of  high  strain  regions,  the  following  extension 
also  applies.  Since  the  exceedance  is  highest  on  Figure  62  it  can  be  assumed 
that  these  regions  necessarily  yielded  first.  The  new  regions  darkened  on 
Figure  60  yielded  intermediately,  and  the  regions  darkened  only  on  Figure  58 
yielded  last.  The  point  is  that  the  second  phase  particles  which  appear  to 
have  connected  on  Figure  58,  have  done  so  by  'growing'  high  strain  regions 
from  the  neighbor  particle,  not  radially  out  from  the  particle  represented  by 
the  center  of  the  particle  model.  Simply  stated,  information  at  the  boundaries 
of  the  particle  model  Is  propogated  inwards  even  though  it  is  harder  to  do, 
than  propogate  plastic  zones  out  from  the  particle,  since  those  elements  tend 
to  be  much  smaller.  It  is  to  be  noted,  that  this  observation  per  se,  does  not 
indicate  any  problem  -  it  is,  however,  interesting. 

In  conclusion,  we  believe  the  model  presented  here,  substantiated  by  the 
results  obtained  and  presented  here,  is  a  tremendous  improvement  over  the  previous 
model.  We  have  paid  great  attention  to  eradicating  the  flaws  of  the  previous 
model,  and  have  succeeded.  We  believe  that  the  model  will  allow  meaningful 
analyses  to  be  performed  that  will  readdress  the  issue  of  interaction  of  second 
phase  particles,  and  its  importance. 
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Figure  3.  Bipolar  element  model,  showing  3  of  5  available  rotations 

of  particle-particle  axis.  Excitation  is  vertical,  constraint 
horizontal. 


Figure  4.  Enlargement  of  central  region 
of  bipolar  map. 
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Figure  6. (for  coherent  model): 

Octahedral  strain  at  inclusion-matrix  interface  as  a  function  of 
angular  position,  four'levels  of  excitation.  Full  constraint, 
rotation  is  0.0  deg.. 
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Figure  7.  (for  coherent  model): 

Octahedral  strain  at  inclusion-matrix  interface  as  a  function  of 
angular  position,  four  levels  of  excitation.  Full  constraint, 
rotation  is  4S.0  deg. 
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Figure  8.  Regions  with  greatest  propensity  for  extensive  strain, 
full  constraint,  rotation  is  0.0  deg,  coherent  model 


.gure  9.  Regions  with  greatest  propensity  for  extensive  strain. 
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Figure  11 
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Enlargement  of  Particle  Model  and  regions 
directly  contguous. 
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1.  ABSTRACT 

This  report  describes  a  simple  substructuring  technique  for  the  finite  element 
analysis  of  elastic  bodies,  one  that  can  be  implemented  by  an  analyst  without  any 
internal  modification  to  an  existing  code.  In  common  with  other  substructuring 
approaches  the  technique  offers  greater  resolution  for  less  computational  effort  than 
would  otherwise  be  required,  and  this  is  demonstrated  by  the  analysis  of  a  test 
problem  involving  a  crack. 

2.  INTRODUCTION 

In  general,  substructuring  for  static  structural  analysis  using  the  finite  element 
method  (FEM)  is  based  on  subdividing  a  large  structure  into  smaller  parts  which  are 
analysed  separately  to  obtain  relationships  between  forces  and  displacements  at  the 
substructure  interfaces.  These  interface  variables  are  then  determined  and  the  results 
used  to  obtain  the  unknowns  within  each  substructure.  Essentially  these  approaches 
involve  partitioning  of  the  stiffness  matrix  in  some  sensible  way  (see  Noor,  Kamel 
and  Fulton  [1]  for  a  fuller  description). 

Our  approach  to  substructuring  consists  of  analysing  an  entire  body  for  stresses 
and  displacements  using  FEM,  then  taking  the  displacements  from  this  first  analysis 
at  a  set  of  interior  nodes  as  the  boundary  conditions  for  a  second  analysis  which 
now  occupies  a  smaller  region.  This  smaller  region  is  centered  on  the  aspect  of  the 
original  problem  of  greatest  interest  and  has  boundaries  which  are  geometrically 
similar  to  those  of  the  original  problem  so  that  the  same  finite  element  map  can  be 
used  in  the  second  analysis.  This  focusing  process  can  then  be  repeated,  with  a 
third  region  within  the  second  being  analyzed  and  so  on*.1  . 

Both  approaches  offer,  in  effect,  the  opportunity  for  greater  resolution  for  the 
same  computational  effort.  This  capability  is  particularly  advantageous  in  analyses 
where  there  are  some  critical  local  regions  such  as  a  crack-tip.  The  main  advantage 
of  the  approach  presented  here  over  those  reported  in,  for  example,  Noor  et.al.  [1], 
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lies  In  its  simplicity:  it  can  be  readily  implemented  by  an  analyst  without 
modification  to  an  existing  coao,  i.e..  externally  to  a  standard  finite  element  program. 

In  what  follows  we  analyze  a  test  problem  for  a  crack  geometry  having  an  exact 
solution  and  show  the  level  of  resolution  that  can  be  achieved  with  a  modest  amount 
of  effort  when  our  substructuring  technique  is  used. 


3.  TEST  PROBLEM 

In  this  section  we  describe  an  elastic  test  problem  with  an  exact  solution  for  a 
configuration  containing  a  crack. 

To  construct  a  test  problem,  we  seek  some  set  of  the  field  quantities-  stresses  <rx, 
try.  tx y  and  displacements  u,  v  -  throughout  a  region  R  satisfying  the  field  equations 
of  two-dimensional  elasticity.  Moreover,  to  provide  a  stringent  test  of  the 
substructuring  procedure,  we  want  these  field  quantities  to  include  a  singularity.  One 
simple  way,  then,  of  constructing  a  suitable  test  problem  is  to  take  the  first  singular 
symmetric  eigenfunction  for  the  crack  as  identified  by  Westergaard  [2],  viz. 
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and  then  to  pick  R  so  as  to  include  the  crack-tip,  and  take  as  boundary  conditions 
on  R  the  values  implied  by  (1).  In  (1),  K  is  a  constant,  G  the  shear  modulus,  v 

Poisson's  retie.  end  r  m'd  9  are  coordi. irtes. 


Fig.  1  Geometry  and  coordinates  for  the  legion  R. 
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A  suitable  R  for  analysis  is  shown  in  Fig.  1,  and  consists  of  a  semi-hexagon  within 
the  upper  half-plate,  tnat  is. 


R  *  {(*»>') | 


-  a  + 


-X  <  x  <  a 


0  <  y  < 


(2) 


where  a  is  the  length  of  a  side  of  the  hexagon.  The  boundary  conditions  on  R  are 
taken  to  be  the  stress-free  conditions  on  the  crack-flank,  the  symmetry  conditions 
ahead  of  the  crack,  and  the  displacements  on  the  remainder  of  the  boundary  as 
prescribed  in  effect  by  (1),  i.e. 

<ry  *  rxy  *  0  on  y  *  0(-a<x<  0), 
v  «  0.  rxy  *  0  on  y  =  0  (0<x<a), 

u  =  J ^  sin  y  cosec  (8  +  -j)  cos  j  [1  -  2v  +  sin2  j]  on  y  =  (a-x) 

v  *  sin  y  cosec  (8  +  y)  sin  y  [2  -  2v  -  cos2  |-]  (0  <  0  <_  ir/3) , 

u  ~  sin  j  cosec  6  cos  j  [1  -  2v  +  sin2  y]  on  y  =  fia/2 
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It  follows  automatically  that  (1)  is  the  exact  solution  to  the  elastic  analysis  of  R 
with  boundary  condition  as  in  (3)  and  thereby  constitutes  a  suitable  test  problem. 


4.  ANALYSIS 

Here  we  describe  the  c.'.c.vS.s  of  the  test  problem  using  our  substructuring 

technique  and  a  constant-strain-triangle,  finite  element  code. 

Since  the  issue  of  interest  is  the  effectiveness  of  our  substructuring  technique,  we 
discretize  R  uniformly  with  a  set  of  equivalent  elements  of  equal  area.  To  examine 
convergence  we  use  a  sequence  of  three  maps  -  a  coarse,  a  medium  and  a  fine  map 
-  with  the  medium  being  constructed  from  the  coarse  by  halving  the  element  size  and 
the  fine  being  similarly  derived  from  the  medium.  The  coarse  map  has  48  elements 
(35  nodes),  the  medium  192  elements  (117  nodes),  the  fine  768  elements  (425  nodes), 
(Fig.  2). 

The  first  step  is  to  analyze  the  original  problem  using  the  sequence  of  maps.  The 
particular  sequence  of  maps  chosen  enables  a  straightforward  but  effective 
extrapolation  technique  to  be  used  to  improve  the  results.  This  extrapolation 
technique  is  based  on  the  following  argument.  We  define  as  our  error  measure  in 
the  displacement  the  difference  e  between  the  values  obtained  by  the  FEM  analyses 
and  the  exact  answer  and  model  this  error  distribution  with 

e  *  e°  h4  ,  (4) 

where  e°  is  the  error  with  a  mesh  size  h  of  unity  and  a  is  a  parameter  reflecting 
the  convergence  of  the  method  (a>0  for  convergence).  For  example,  the  errors  in  the 
displacements  in  the  x-direction  are 


c  c  »  o,  a 
e  =  u-  u  =  eh. 


m  m  «  o  .h, 

e  =  u  -  u  =  e  , 

f  f  00  o..h..a 

e  *  u  -  u  =  e  (^)  , 


where  the  superscripts  c,m,f  and  »  denotes  respectively  the  coarse,  medium, 
fine,  and  "exact"  (or  infinite  number  of  elements)  results.  Eliminating  e°  and  a  from 
(5)  then  gives 


5MW 
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similarly. 


(7) 


Equations  (6)  and  (7)  are  extrapolated  estimates  of  u  and  v  which  realize 
improvements  over  the  fine  map  values. 

The  next  step  is  to  perform  a  first  substructure.  We  select  as  our  subregion  R,  a 
geometrically  similar  region  to  R  centered  on  the  crack-tip  with  sides  equal  to  half 
those  of  R:  thus 

R*  «*.y)|  -!♦-*<*<!--*,  0<  v<%- 

rHS  /3 

(8) 

To  obtain  the  boundary  conditions  on  "R  from  the  analysis  of  R,  we  proceed  as 
follows.  We  use  the  estimates  u”  ^obtained  from  (6)  and  (7)  at  nodes  where  there  are 
ue,  um.  uf  and  vc,  vm.  vf  values.  Where  u®  and  v”  can  not  be  found  in  this  manner 
because  uc  and  vc  are  not  available,  we  interpolate  between  the  u®  and 

CP  00 

v  determined  from  (6)  and  (7)  using  the  fine  map  distribution.  That  is,  let  u}  ,  and 
u5  be  the  extrapolated  values  of  u  at  nodes  1  and  5  and  suppose  that  the  nodes 
of  the  intervening  three  fine  map  values  are  not  supported  by  coarse  map  results. 

CD  00  ft 

We  estimate  these  intermediate  node  values  --  u  2  *  u^  ,  ,  by  fitting  the 

corresponding  distribution  of  fine  map  values  to  the  end  points  u^  ,  u*  in 
accordance  with2  . 
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where  ai  and  a  2  are  defined  as 


(10) 


In  instances  where  fine  estimates  are  not  even  available  -  as  at  some  ndoes  in  the 
fine  map  for  TT-  we  linearly  interpolate  over  U.“  and  V."  determined  as  described 
above. 

Using  these  extrapolated  and  interpolated  values  as  boundary  conditions  and 
analyzing  R  constitutes  our  first  substructure:  applying  the  same  procedure  twice  is 
the  second  and  so  on.  In  actuality,  only  two  substructures  were  performed  using  the 
constant-  strain-triangle,  elastic,  finite  element  code  PLANDJ  [3].  That  is,  we  ran 
the  coarse  medium  and  fine  maps  for  the  pre-substructure,  first  substructure  and 
second  sub-structure  geometries.  In  this  manner  we  hoped  to  enjoy  the  resolution  of 

3 

"effective"  coarse,  medium  and  fine  maps  having  120  elements  (48  *  x  *  48  ♦ 

3  H 

•  48},  480  elements  and  1920  elements,  respectively  (Fig.3  sketches  the  effective 

fine  map).  Our  storage  requirements  using  substructuring  however  are  about  40%  of 

those  required  for  a  direct  analysis  of  these  effective  maps  and  the  number  of 

operations  needed  using  substructuring  is  approximately  1/3  of  the  number  required 

for  a  direct  analysis  even  when  the  multiple  analyses  required  for  the  substructuring 
are  taken  into  account.3 
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where  N  is  the  number  of  degrees  ox  freedom  and  B  is  the  band-width 
Zienkiewicz  [4],  p.463. 
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S.  RESULTS 

In  this  section  -.ve  show  tha  results  from  the  application  of  our  substructurinr 

technique,  and  in  particular  assess  the  increased  resolution  gained  at  the  crack-tip 
using  substructuring. 

Figs.  4  through  6  show  the  convergence  of  the  initial  analysis  by  comparing  the 

coarse,  medium  and  fine  results  with  the  exact  solution  of  (1).  In  Figs.  4  and  5  the 

dimensionless  stresses,  a  -  /a"o  IK  and  a  »  ^ o  IK  formed  from  nodal  averages, 

xx  y  i_ 

are  plotted  as  a  function  of  dimensionless  distance  r  =  r/a  along  a  ray  at  8  «  W 3. 
The  stresses  converge  except  for  the  peak  value  at  the  innermost  node,  i.e.  nearest 
the  crack-tip  (rsO).  Fig.  6  shows  the  dimensionless,  crack-opening,  nodal  displacement, 
v*Gv/K  fi  at  0  *  »,  as  a  function  of  T.  The  displacements  converge  even  at  the 
crack-tip.  This  convergence  allows  us  to  use  the  extrapolation  technique  outlined 
earlier  to  form  the  boundary  conditions  in  our  first  substructure. 

Figs.7  through  10  show  the  displacement  values  used  as  the  boundary  conditions  in 
the  application  of  the  substructuring  technique,  comparing  the  coarse,  medium  and 
fine  map  results  with  the  exact  solution  of  (1).  In  Figs.  7  and  8  the  dimensionless 
nodal  displacements,  u  *  Gu/K  ^a  and  ~v~  *  Gv/K  ^a,  around  the  upper  boundary  of  R 
are  shown  i.e.  as  functions  of  the  dimensionless  distance  s  *  s/a  from  the  crack  8  = 
*  to  8  •  0.  while  the  corresponding  displacements  obtained  from  the  analysis  of  R 
for  use  in  the  second  substructure  are  shown  in  Figs.  9  and  10.  For  both  sets  of 
curves  the  displacement  values  derived  using  the  extrapolation  interpolation  technique 
described  in  the  previous  section  are  not  distinquishable  from  the  exact  values  on  the 
size  scale  used  in  the  figures.  Accordingly  no  accumulation  of  error  in  the  boundary 
values  results  from  the  present  multiple  application  of  our  substructuring  technique. 

Lastly,  we  examine  the  error  in  the  crack  opening  displacement  near  the  crack-tip 
as  the  substructuring  proceeds.  Since  v  at  8  *  n  goes  to  zero  at  7^0  we  normalize 
this  error  by  dividing  by  a  typical  specific  value  of  v(0*r)  rather  than  the  exact 
crack-opening  value  at  the  corresponding  7  -  position:  the  actual  typical  value  taken 
is  the  crack-opening  displacement  at  the  crack  mouth  in  R,  that  is  at  r  *  1.  The 
error  measure  in  extrapolated/interpolated  values  is  then  given  by 

e  •  (v-  -  v*)/v* 

where  v  is  the  displacement  found  by  extrapolating/interpolating  the  FEM  results 
at  8*n.  as  distinct  from  the  exact  soluticr  ve  there,  and  vl  in  the  typical  n? 
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Fig.  9  Displacement  (u)  along  the  boundary  after  first  substructuring 


Displacement  (v)  along  the  boundary  after  first  substructuring 


v.  Fig.  11  shows  e  after  the  first  analysis  (pre-substructuring)  and  after  the  first  and 
second  substructure.  Tne  peat,  error  occurs  at  the  node  adjacent  to  tne  cracr.  tip  m 
all  three  analyses  and  is  uniformly  reduced  with  substructuring.  An  alternative  way 
of  viewing  this  improvement  is  to  track  the  reduction  in  extent  of  the  region  near 
the  crack-tip  in  excess  of  some  error.  Taking,  for  example,  the  region  in  which  e  £ 
0.01  then  the  associated  radial  dimensions  are  0.282,  0.107  and  0.037  from  the  pre-, 
first,  and  second  substructure  analyses  demonstrating  the  added  resolution  gained  by 
the  technique  near  the  crack-tip.  In  practice,  too,  if  it  was  necessary  to  drive  this 
region  to  zero  at  this  error  level,  further  substructures  could  be  performed:  thus  the 
technique  can  potentially  be  adjusted  to  achieve  whatever  level  of  resolution  is 
desired. 


Pre - subs t  ructur in g 
First  substructuring 

Second  substructuring 


Fig.  11  Error  in  the  crack  opening  displacement. 
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COMPUTATION  OF  SLOW,  STABLE  CRACK  GROWTH 
IN  A  FRACTURE  TEST  SPECIMEN 
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Pittsburgh,  Pennsylvania 

and 

C.W.  Cho 
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ABSTRACT 

The  title  problem  has  been  solved  using  a  compact  specimen  of 
7475-T6  aluminum.  It  is  found  that  loading  forces  are  accurately  predicted. 
Also,  strain  at  the  crack  tip  undergoes  a  modest  rise  during  the  course 
of  crack  advance  so  that  constant  crack  tip  strain  is  inconsistent  with 
J-controlled  crack  growth.  That  these  results  are  consistent  with 
experimental  measurements  leads  to  some  discussion  of  a  number  of  computational 
sound  round  robins  recently  reported  in  the  literature. 


V' 


FILMED 


